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1.0  INTRODUCTION 


While  a  wide  variety  of  pn^ulsion  systems  have  been  developed  over  the  years,  a 
gap  has  developed  in  the  medium  range  of  specific  impulse  (roughly  1(XX)  to  2000  s)  for 
which  no  advanced  electric-propulsion  system  is  ready.  Chemical  rockets  cover  the 
territtny  below  3(X)  s  effectively,  and  state-of-the-art  ion  thrusters  are  highly  efficient  in  the 
2:3(XX>-s  regime. 

The  ccmcept  of  the  orbit-transfer  vehicle  (OTV)  has  evolved  in  recent  times;  these 
vehicles,  which  would  carry  spacecraft  from  low  orbits  to  higher  orbits,  would  be 
t^timally  served  by  a  propulsion  system  with  a  specific  impulse  Isp  in  the  16(X)-s  range. 
(Appendix  1  lists  all  symbols  used  in  this  report)  The  OTV  peifonnance  requirement  that 
drives  this  optimum  is  the  need  to  achieve  short  trip  times  with  the  available  onboard 
electrical  power,  that  is,  the  ratio  thrust/)x>wa'  F/P  is  determined  by  the  available  power  P 
and  the  thrust  F  needed  to  achieve  the  desired  trip  dme.  This  ratio  of  thrust  to  power  must 
be  provided  with  minimum  propulsion-system  mass,  to  avoid  degrading  the  payload- 
maWspacecraft-mass  ratio.  Ihe  two  key  relatitmships  are 

F/p  =  ^  (1.0-1) 

and 

(1.0-2) 

where  1)7  is  die  total  efficiency  of  the  thrust  system,  g  is  the  acceleration  due  to  gravity.  It 
is  the  total  mission  impulse  requirement  (total  dirust  x  dme),  and  Mp  is  the  total  prc^llant 
mass.  These  relationships  show  the  importance  of  specific  impulse  in  determining 
propellant  mass  and  specific  impulse  combined  with  overall  thruster  efficiency  in 
determining  dunst-to-power  ratio.  The  equations  show  the  conflicting  trends  toward  high 
specific  impulse  to  minimue  propellant  mass  and  low  specific  impulse  to  maximize  thrust- 
to-power  ratio  and  therefore  to  minimize  trip  time. 

Equatxxis  1.01-1  and  1.01-2  can  be  combined  to  find 
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(1.0-3) 


that  is,  for  a  given  thrust-to-power  requirement,  the  total  propellant  mass  is  inversely 
proportional  to  the  total  thrust-system  efficiency,  qj;  if  the  total  efficiency  is  50%,  then 
twice  the  propellant  mass  is  needed,  compared  to  an  ideal  thruster. 

Several  novel  thrust  systems  have  been  developed  in  response  to  the  OTV 
requirement  for  medium  specific  impulse.  Rocket  Research  Company,  for  example,  has 
developed  a  flightready  arcjet  thruster  that  will  be  flown  on  the  G£  Astrospace  Series  70(X) 
satellite,  beginning  in  1993.  The  6-kg-mass  thruster  operates  on  a  power  of  1.8  kW  to 
produce  thrust  at  specific  impulse  of  ^00  s.  A  similar  but  larger,  26-kW  ammonia  arcjet 
operates  at  a  specific  impulse  of  812  s,  and  is  expected  to  be  flight  tested  in  1995.  llie 
former  Soviet  Union  had  developed  a  “stationary-plasma”  thruster  operating  at  1600  s  of 
specific  impulse,  providing  up  to  350  mN  of  thrust  at  a  6-kW  power  level.  Still,  these 
devices  have  certainly  not  closed  the  specific-impulse  gap:  the  arcjets  still  operate  at  an 
undesirably  low  level  of  specific  impulse,  and  the  Russian  device  appears  to  experience 
severe  life-limiting  erosion  processes. 

Conventional  ion  thrusters,  employing  inert-gas  propellants  such  as  xenon,  can 
operate  at  low  specific  impulse  simply  by  using  reduced  beam  voltages.  However,  two 
problems  arise  which  have  thus  far  prevented  the  development  of  capable  thrusters  in  the 
1000-to-2(X)0-s  specific-impulse  range:  first,  the  thruster  electrical  efficiency  tie  falls  off 
with  low  extraction  voltages;  and  second,  operation  at  low  beam  voltage  requires  high 
beam  currents  to  provide  unreduced  thrust,  but  the  high  beam  currents  cause  early  wearout 
of  the  extraction  grids.  If  a  heavy  propellant  particle  were  available,  thrust  could  be 
provided  at  higher  beam  voltages  and  Iowa*  beam  currents,  and  these  problems  would  be 
averted.  Ousters,  with  their  high  mass  compared  to  monoatomic  inert  gases,  offer  the 
promise  of  filling  this  need. 

Figure  1  shows  our  concept  for  ion  propulsion  using  ionized  clusters.  The  basic 
idea  is  that  atomic  or  molecular  species  introduced  into  the  nozzle  shown  on  the  left  will, 
under  conditions  discussed  in  the  next  section,  condense  into  clusters  ranging  from  a  few 
hundred  to  a  few  thousand  atoms.  These  clusters  can  be  readily  ionized  and  accelerated  to 
high  energy,  just  as  if  they  were  heavy  irxrnoatomic  inen  gases. 
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Figure  1.  Schematic  of  the  ionized-cluster  apparatus  used  for  this  smdy. 

Under  this  contract,  we  have  made  use  of  an  apparatus  constructed  under  a 
previous  effort^  to  test  the  idea  of  using  clusters  to  form  a  heavy  propellant  ion  for  space- 
propulsion  qjplications.  Basically,  we  have  been  successful  in  this  attempt,  although  a 
great  deal  of  additional  work  remains  to  be  done  to  bring  the  concept  to  application.  We 
have  been  able  to  demonstrate  thrust  levels  above  10  mN  in  our  relatively  primitive 
labcMiattny  q)paratus,  and  we  have  achieved  near-oxnplete  clusterization  of  carbon-dioxide 
gas,  when  supplied  with  hydrogen  carrier  gas  to  aid  in  the  cluster-condensation  process. 
Our  work  buil^  upon  earlier  studies^  in  achieving  high  clusterization  with  molecular  gases; 
molecular  gases  fmm  smaller  clusters  with  masses  around  15,000  amu,  compared  to  inert- 
gas  clusters,  which  typically  have  masses  of  the  order  of  250,000  amu.  The  smaller-mass 
clusters  allow  beam  fmmaticni  with  lower  voltages,  vdiich  may  be  a  practical  advantage. 

The  key  processes  that  appear  to  block  the  patii  of  itxtized-cluster  beams  as  a  means 
of  q)ace  propulsion  are  the  need  for  excessive  amounts  of  carrier  gas  (which  degrades  the 
prc^llant-utilization  efficiency),  the  pom*  performance  of  the  ionizer  (which  requires 
excessive  power  and  produces  excessively  wide  distributions  of  mass-to-charge  ratio)  and 
the  beam-extraction  system,  which  accelerates  cmly  a  small  faction  of  the  incident  ionized- 
cluster  current  We  telieve  that  all  of  these  difficulties  may  yield  to  further  investigation, 
and  the  preliminaty  accomplishments  that  we  have  made  in  this  study  may  lead  to  a  viable 
^pace-prt^ulsi<m  sysrem. 
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2.0  BASIC  CLUSTER-BEAM  PHYSICS 


In  this  section  we  describe  the  basic  processes  by  which  ionized  cluster  beams  are 
formed:  cluster  nucleadon  and  growth  with  supersonic  nozzles  and  carrier  gases,  cluster 
ionization  with  electron  bombardment,  and  the  extraction  and  formation  of  an  ionized- 
cluster  beam.  We  will  also  discuss  the  characterisdcs  of  ion-propulsion  systons  so  that  we 
can  evaluate  the  potential  of  ionized  cluster  beams  for  space-propulsion  applications. 


2.1  CLUSTER  FORMATION 

In  this  section  we  present  a  simple  model  which  illustrates  the  principal  elements  of 
the  phenomena  that  are  observed  in  gas  expansion  through  supersonic  nozzles  under  a 
specific  set  of  conditions.  This  model  is  similar  to  others^*^  that  have  appeared  previously 
in  die  literature. 

Ousters  are  formed  when  gas  expanding  through  a  supersonic  nozzle  cools  m  the 
supersaturation  point;  it  then  begins  to  nucleate  and  condense,  forming  clusters.  Figure  2 
illustrates  the  basic  geometry  for  cluster  formation  in  a  supersonic  nozde.  Gas  is  supplied 
to  the  nozzle  at  the  stagnation  pressure,  which  is  the  pressure  corresponding  to  the 
nonflowing  gas  within  the  supply  reservoir.  If  the  gas  expands  without  heat  addititm  or 
fiictional  losses  (i.e.,  the  expansion  is  isentropic),  then  the  internal  energy  and  pressure  of 
the  gas  is  converted  into  kinetic  motion  of  the  atoms  or  molecules.  During  this  isentropic 
expansion  the  temperature  of  die  gas  decreases,  following  the  curve  shown  in  Fig.  3.  Fot 
id^  gases,  the  isentrope  can  be  characterized  by  the  equation 

TV^‘ =  const.  ,  (2.1-1) 

where  v=Cp/cv.  and  Cp  and  Cy  are  the  speciHc  heats  of  the  gas  at  constant  pressure  and 
volume,  respectively.  When  the  pressure  and  temperature  of  the  expanding  gas  pass  to  the 
left  of  die  vapor-pressure  curve  (refer  to  Fig.  3)  fm*  die  particular  gas  in  questicm,  it  enters  a 
region  of  supersatuiation,  i.e.,  the  vapor  density  is  above  that  which  can  be  sustained  at  die 
given  temperature.  In  this  transient  condition,  the  gas  will  begin  to  condense,  forming 
nuclei  which  eventually  grow  into  clusters.  During  condensation,  however,  the  latent  heat 
of  condensation  is  released,  which  briefly  increases  the  temperature  and  pressure  of  the 
expanding  gas  as  illustrated  in  Fig.  2.  This  increase  in  temperature  brings  the  gas  back  to 
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the  viyx>r-pFessure  curve,  and  acts  as  feedback  mechanism,  preventing  further  excursions 
into  the  supersaturadon  regitm. 

According  to  Reference  (3),  cluster  nuclei  are  formed  during  the  brief  initial 
excursion  into  the  supersaturadon  region.  These  clusters  wiU  condnue  to  grow  as  long  as 
the  collision  frequency  with  unclusterized  atoms  or  molecules  is  significant,  but  addidonal 
clusters  are  not  formed. 
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Figure  2.  Pressure  and  temperature  distribudons  in  gas  flow  through  a  supersonic  nozzle 
with  condensadon.  According  to  Oswadtsch^  nucleadon  takes  place  only  within  a  narrow 
zone,  while  cluster  growth  continues  altmg  much  of  the  nozzle. 
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Figure  3.  Pressure-temperature  diagram  of  gas  expansion  in  a  supersonic  nozzle  with 
condensation.  As  the  gas  enters  the  nozzle  it  exp^s  isentropically  at  first.  When  it 
reaches  the  point  of  saturation,  condensation  can  set  in.  On  its  further  expansion,  the  gas 
tracks  the  pressure  curve  and  absorbs  the  latent  heat  of  condensation  as  it  proceeds. 

2.2  CLUSTER  IONIZATION 

Cluster  ionization  is  obtained  in  our  experimental  apparatus  by  neutral-cluster 
bombardment  with  electrons  derived  from  a  ring-shaped  filamait  and  accderated  throu^  a 
cylindrical  screen-mesh  anode. 

The  electron  current  in  the  ionizer  required  to  assure  near-complete  ionization  of  the 
neutral  cluster  beam  can  be  estimated  from  a  single  collision-cross-section  calculation: 


L  a,N’ 


(2.2-1) 


where  N  •  300  to  600  is  the  number  of  atoms  in  a  cluster,  vd  •  10^  tn/s  is  the  cluster 
velocity,  Oj «  2  to  3  x  lO'^O  are  typical  ionization  cross  sections,  and  using  an  ionizer 
axial  length  of  L  =  20  mm,  we  find  that  current  densities  of  the  ordo-  of  je  •  1  kA/m^ 
(100  mA/cm2)  are  needed  to  ensure  a  high  itmization  fraction.  (Experimentally,  we  have 
found  that  increasing  die  electron  current  density  causes  the  cluster-beam  current  to  increase 
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until  electror  current  densities  near  the  1  kA/m^  value  are  reached;  this  observation  seems 
to  confir.^  <Mir  rough  calculation.) 

In  order  to  extraa  the  magnitude  of  electron  current  just  determined,  we  must 
design  an  ionizer  within  the  constraints  of  the  Qtild-Langmuir  space-charge  limit  on  current 
density: 


where  ZQ  is  the  cathode-to-anode  gap,  eq  is  the  permeability  of  free  space,  Vq  is  the 
electron  accelerating  voltage  (discharge  voltage),  and  e  and  me  are  the  charge  and  mass  of 
the  electmi,  respectively.  Using  a  typical  value  of  Vd  =  100  V,  we  find  that  we  require  a 
gap  ZQ  ^  1.6  pm  in  order  to  achieve  the  required  1  kA/m^.  Maintaining  such  a  small  gap 
between  two  material  electrodes  is  obviously  impractical  (particularly  when  one — the 
filament — must  be  heated  to  roughly  a  2000-K  temperature).  However,  the  presence  of 
ionized  atomic  species  (ions  formed  from  unclustered  thrust  gas  and  carrier  gas  atoms) 
leads  to  the  formation  of  a  plasma  sheath  between  the  anode  and  cathode  which  greatly 
narrows  the  effective  anode-cathode  gap  below  the  physical  separation  of  the  met^ 
electrodes. 

As  is  well  known  in  plasnui  diodes  similar  to  our  ionizer,  the  bulk  of  the  plasma 
remains  near  anode  potential,  and  the  100-V  discharge  voltage  appears  across  a  thin  sheath 
near  the  cathode.  The  characteristic  thickness  of  the  sheath  is  the  plasma  Debye  length, 
which  is  in  the  40-pm  range  in  plasmas  of  modest  density  (ne  «  10^^  m*^).  We  believe  that 
the  presence  of  this  ionizer-filling  plasma,  formed  of  the  waste  (non-cluster)  gases 
emanating  from  the  supersonic  nozzle,  serves  two  important  functions:  it  permits  extractitm 
of  high  electron-current  densities,  as  just  described,  and  it  also  causes  nearly  all  of  the 
clusters  to  become  ionized  at  the  same  electrostatic  potential — that  of  the  anode  plasma. 

The  foregoing  observations  pdnt  to  the  need  for  further  work  in  two  specific  areas: 
improved  ionizer  design  to  exploit  plasma  effects,  and  adjustment  of  the  electron  energy 
eVo  to  provide  as  uniform  as  possible  a  mass-to-charge  ratio,  Mci/q.  In  die  case  of  ionizer 
design,  we  feel  that  incorporating  modem  magnetic-confinement  methods  (which  have 
been  highly  developed  in  Penning-discharge  ion  thrusters)  could  dramatically  improve 
kmuer  efficiency  in  tenm  of  the  ionizer  power  required  to  ionize  an  an:^)ere  of  cluster-beam 
current  In  die  case  of  the  mass-to-charge  ratio,  we  will  show  in  a  later  section  that  having 
a  broad  qnead  in  MciAl  significandy  degrades  the  efficiency  of  any  cluster-beam  thruster 
that  could  be  designed.  While  sotiK  ^nead  in  Mci/q  is  inevitable  due  to  the  randomness  of 
the  clusterization  [ffocess,  die  fact  that  the  ionizer  creates  clusto’  ions  with  charges  ranging 
from  q  s  e  to  4e  imposes  an  unacceptable  loss  in  thrusting  efficiency.  While  we  have  not 
yet  atusn^ited  to  reduce  the  range  of  cluster  charges,  we  feel  that  a  narrower  distribution 
could  be  achieved  by  combining  reduced  electron  energy  (around  40  eV)  with  the  greater 
ionizer  control  and  ^ciency  afforded  by  the  Penning  discharge. 
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2.3  CLUSTER-BEAM  EXTRACTION  SYSTEM 


The  cluster-beam  extraction  system  that  we  have  employed  is  a  derivative  of 
electrostatic  beam-forming  technology  that  has  been  developed  for  application  in  ion 
thrusters  and  broad-beam  industrial  ion  sources  over  the  last  30  years,  with  a  simple 
modification  to  exclude  efficiency-degrading  light  ions  ^m  the  beam.  We  want  our 
extraction  system  to  form  its  beam  without  extracting  light  ions  and  to  accelerate  heavier 
ions  without  signiHcant  interception  of  ion  trajectories  by  the  grids  themselves.  When  there 
is  interception,  we  require  that  the  energy  per  incident  ion  be  below  the  sputter  threshold^ 
so  that  the  grids  will  not  be  rapidly  eroded  away. 

A  conunon  embodiment  of  the  broad-beam  extraction  system  employs  three  grids, 
each  fabricated  from  a  thin  refractory  metal  (commonly,  molybdenum)  with  a  densely 
packed  hexagonal  pattern  of  apertures.  The  apertures  between  the  three  grids  are  aligned  in 
such  a  manner  that  the  emerging  beamlets  are  parallel  to  the  thrust  axis.  In  the  case  of  flat 
grids,  the  apertures  are  paraxial;  if  the  grids  are  formed  into  spherical  caps  for  mechanical 
stability  (as  is  common  practice),  then  the  apertures  must  be  slightly  displaced  or 
compensated  to  provide  paraxial  beamlets  and  avoid  cosine-law  thrust  loss.  The  three  grids 
are  biased  as  follows:  the  innermost  or  screen  grid  is  maintained  at  the  high  positive  beam 
voltage,  the  intermediate  or  accel  grid  is  held  negative  (typically  a  few  hundred  volts)  to 
prevent  electrons  from  outside  the  thruster  from  flowing  upstream  into  the  thruster,  and  the 
outermost  or  decel  grid  is  at  ground  potential. 

Ion  trajectories  in  a  three-grid  beam  extraction  system  vary  primarily  because  of 
differences  in  the  point  of  entry  into  the  beam-extraction  system  and  because  of  ion- 
temperature  effects.  In  our  cluster  beams,  however,  the  mass-to-charge  ratio  varies 
widely,  due  to  the  formation  of  multiply  charged  clusters  and  to  the  statistical  variation  in 
cluster  masses  formed  in  the  supersonic  nozzle.  This  mass  variation,  which  has  been 
observed  to  be  gaussian  in  character,^  will  potentially  lower  the  thrust  density  because  the 
species  having  higher  mass-to-charge  ratio  will  cause  the  extraction  to  become  perveance 
limited  prematurely.  Another  potentially  troubling  problem  is  that  many  carrier-gas  atoms 
become  ionized;  their  masses  of  2  to  4  amu  compare  with  20  to  50  x  10^  amu  for  molecular 
and  inert-gas  clusters.  These  light  ions  pose  a  twofold  problem:  they  seriously  degrade 
thruster  performance,  as  noted  above,  and  they  can  cause  extraction- system  design  to 
become  ^fficult  because  of  the  very  different  initial  velocities  of  the  cluster  ions  and  the 
light  ions.  Fortunately,  we  have  a  simple  solution  to  the  light-ion  problem:  our  ion- 
extraction  system  incorpOTates  a  fourth  grid  inside  the  screen  grid  to  reject  light  ions.  If  the 
potential  on  the  retarding  grid  is  Vr,  then  clusters  having  a  mass-to-charge  ratio  of 


OT  mwc  will  be  accepted  into  the  extraction  system;  lighter  ions  or  clusters  will  be  rejected. 
The  comparatively  small  fractitm  of  these  lighter  ions  makes  the  mass-utilization  cost  of  this 
apjntMch  a  small  one. 
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2.3.1  Extraction-System  Sputtering  by  Cluster  Ions 

The  remaining  beam-extraction  issue  confronting  a  practical  ionized-cluster  thruster 
is  sputtering  of  extraction  grids  by  clusters  themselves.  First  of  all,  there  are  three  types  of 
sputtering  to  consider,  corresponding  to  the  specific  grid  involved : 

Type  1.  Screen-Grid  Sputtering  by  impact  of  those  clusters  that  strike  the 
upstream  side  of  the  screen  grid. 

Type  2.  Charge-Exchange  Sputtering  of  the  accel  grid  by  ions  that  are  formed 
when  ionized  clusters  strike  neutral  gas  atoms  (carrier  gas  or 
unclusterized  thrust  gas)  and  transfer  some  charge  to  the  neutral,  and 

Type  3.  Direct  Interception  of  ionized  clusters  that  strike  grids  after  they  have 
been  accelerated  to  high  energy. 

Sputtering  of  the  fourth  (retarding)  grid  is  unlikely  to  occur,  since  it  is  biased 
positive  with  respect  to  anode  voltage,  and  the  decel  is  unlikely  to  be  sputtered  by  charge- 
exchange  ions  because  of  their  low  energy.  We  discuss  these  three  types  separately  below. 

Type-1  sputtering  is  trivially  limited  by  keeping  the  energy  of  the  incident 
ionized  clusters  below  the  sputter  threshold  of  the  screen  grid.  Since  most  of  this  energy  is 
acquired  through  acceleration  from  the  ionizer-plasma  potential  to  the  screen,  it  is  trivially 
controlled  by  maintaining  ionizer  plasma  potentials  low.  In  our  experiments,  we  have  no 
ability  to  investigate  this  approach,  since  our  ionizer  does  not  employ  any  form  of  magnetic 
confinemoit  and  achieves  efficient  ionization  only  when  the  voltages  are  near  1(X)  V. 

Type-2  sputtering  normally  results  from  the  fact  that  the  accel  grid  must  be  held 
negative  several  hundred  volts  to  prevent  entiy  of  elections  into  the  extracticm  system  from 
outside  the  thruster.  (Tharge-exchange  ions  formed  within  the  beam-extraction  system  and 
outside  the  thruster  fall  through  all  or  a  part  of  this  potential  difference,  causing  sputter 
erosion  of  accel  surfaces  and  the  walls  of  the  extraction  apertures.  In  conventional  ion 
propulsion,  accel  erosion  can  be  reduced  but  not  eliminated  by  the  use  of  a  outer  decel 
electrode,  which  partially  shields  charge-exchange  ions  from  the  E-fields  associated  with 
the  accel  electrode.  Accel  sputtering  is  a  life-limiting  process  which,  so  far,  has  only  been 
partially  eliminated.  In  ionized-cluster  thrusters,  we  expect  that  accel  erosion  may  be  an 
even  more  serious  problem  than  it  is  in  ion  thrusters,  b^use  of  the  presence  of  the  very 
large  amount  of  carrier  gas.  This  large  amount  of  carrier  gas  will  lead  to  large  charge- 
exchange  currents,  but  the  rate  of  sputtering  will  also  be  reduced  somewhat  because  the 
carrier-gas  ions  (H2  or  He)  are  light  and  have  low  sputtering  yields. 

It  is  also  possible  that  the  wide  spread  in  charge-to-mass  ratios  present  in  the  cluster 
beam  may  associate  with  a  large  spread  in  initial  cluster  velocities.  If  true,  this  could 
require  the  use  of  larger  decel  apertures  to  accommodate  the  varying  particle  trajectories; 
lai^  decel  apertures  would  exacerbate  the  problem  of  accel  erosion  by  allowing  accel 
fields  to  penetrate  into  the  charge-exchange  plasma  downstream  of  the  decel  grid.  Only 
when  further  woiit  has  been  done  will  we  be  able  to  determine  the  magnitude  of  this 
problem  and  design  appropriate  means  to  control  it 
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Type>3  sputtering,  or  direct  interception  on  the  accel  grid  by  clusters  with  the 
energy  imparted  by  the  total  screen-to-accel  voltage  is  caused  by  two  processes;  the  Hrst 
results  from  the  fact  that  clusters  entering  the  extraction  system  have  differing  energies, 
masses,  and  entry  points  into  the  extraction  system.  The  second  results  from  the  fact  that 
the  detailed  shape  of  the  electric  fields  near  the  entry  point  to  the  extraction  system  are 
subject  to  fluctuations  induced  by  collective  plasma  effects.  This  means  that  even  if  all 
possible  cluster  trajectories  will  traverse  the  extraction  system  without  impact  during 
quiescent  operation,  plasma  fluctuations  may  lead  to  new  trajectories  that  will  impact  the 
grids.  This  sputtering  process  would  likely  be  no  more  of  a  problem  than  it  is  in 
conventional  ion  thrusters,  except  for  the  perturbing  effects  of  the  retarding  grid  on  initial 
cluster  trajecmries.  The  next  few  paragraphs,  which  identify  a  maximum  specific  impulse 
that  can  be  achieved  if  we  are  absolutely  to  hold  cluster  energy  below  the  sputter  threshold, 
should  be  viewed  as  speculative,  since  no  such  limitation  is  needed  in  conventional  ion 
thrusters. 

In  conventional  ion  thrusters,  both  of  these  types  of  sputtering  can  be  controlled  or 
eliminated  by  adjusting  extraction-system  design  parameters  and  potentials  within  the 
discharge  plasma.  With  a  cluster  beam,  however,  the  wide  range  of  initial  cluster  velocities 
may  present  difficulties  in  employing  the  conventional  technology.  In  our  studies  we  have 
not  yet  addressed  these  extraction-system  design  issues,  but  we  do  have  a  simple  approach 
to  eliminating  the  direct-interception:  limit  the  cluster  energies  such  that  individual  members 
of  the  cluster  have  energies  below  the  sputtering  threshold. 

The  problem  of  sputtering  of  extraction  grids  by  clusters  has  been  investigated 
previously,^  with  the  finding  that  sputtering  does  not  appear  to  occur  if  the  energy  Wa  of  an 
individutd  atom  within  the  cluster  (viz.,  the  kinetic  energy  that  it  would  have  if  it  had  the 
same  mass  and  velocity,  but  were  not  a  part  of  a  cluster)  is  below  the  sputter  threshold  for 
bombardment  of  the  specific  grid  material  by  the  cluster  atom.  That  is,  if,  for  example,  30- 
eV-per-atom  xenon  atoms  or  ions  will  not  sputter  molybdenum,  then  reference  (S)  suggests 
that  neither  will  any  number  of  them  when  they  are  packed  together  into  a  cluster.  Since 
the  atomic  energy  Wa  can  be  written  as 


W. 


IP 

2R 


(2.3-2) 


the  sputter  threshold  Walthresh  sets  a  limit  on  achievable  sputter-free  specific  impulse, 
independent  of  the  size  of  the  clusters  themselves.  For  xenon,  the  sputter  threshold 
Walthresh  on  molybdenum  is  about  35  eV,  leading  to  a  maximum  specific  impulse  Ispimax 
without  sputtering  of  584  s  (where  we  have  somewhat  arbitrarily  taken  i1t=1  as  a  total 
efficiency  and  R=Vb/Vt=0.8  as  an  upper  bound  on  the  net-to-totsd  voltage  ratio).  This  is 
the  maximum  achievable  sputter-free  speciflc  impulse,  which  may  yet  be  degraded  by 
factors  such  as  beam  divergence,  nonunifoim  mass/charge,  etc. 

Equation  2.3-2  can  be  multiplied  by  the  number  of  atoms  in  a  cluster,  n,  to  yield 
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(2.3-3) 


where  nm  =  Mci  is  the  total  cluster  mass.  Figure  4  shows  cluster-atom  energy  given  atomic 
mass  and  maximum  specific  impulse.  By  choosing  a  sputtering  threshold  voltage  and  a 
value  of  specific  impulse,  one  can  estimate  the  range  of  atomic  masses  that  can  be  used 
without  risking  sputter  damage  to  the  extraction  system. 


Figure  4.  Kinetic  energies  of  atoms  or  molecules  in  clusters  as  a  function  of  the  atomic 
mass  number  for  different  specific-impulse  values.  To  avoid  sputtering  of  the  accel 
electrode,  the  energy  pa  atom  should  remain  below  die  sputter  threshold  (typicaUy  25  to  35 
eV),  and  the  atomic  ituss  number  should  be  limited  accmdingly.  The  vertical  dotted  lines 
conesptmd  (in  increasing  mass)  to  the  masses  of  methane,  ethylene,  argon,  and  carbon 
dioxitte,  respectively. 


2.3.2  Relationship  Between  Cluster  Mass  and  Beam  Voltage 
Equation  2.3-3  can  be  rewritten  in  terms  of  the  beam  voltage  Vg: 

v,=^(l,g)*  (2.3-4) 
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Equation  2.3-4  reveals  that  for  a  given  specific  inqiulse,  the  beam  voltage  Vb  is 
prc^xntional  to  the  cluster  mass.  Figure  5  illustrates  this  rel^onship,  and  shows  the  effect 
of  making  an  arbitraiy  choice  of  upper  bound  of  desired  beam  voltages.  While  very  high 
voltages  can  certainly  be  accomiiKxlated  in  space,  questions  of  power-supply  efficiency, 
weight,  and  difficulty  of  insulaticm  become  significant  as  higher  beam  voltages  are  used. 


Figure  5.  Quster-beam  voltage  as  a  function  of  cluster  mass  for  various  values  of  specific 
inqiulse.  A  particular  (arbitrary)  choice  for  the  upper  bound  of  reasonable  beam  voltages  is 
shown  to  iUustrate  the  effect  on  choosing  an  apivopriate  cluster  mass. 


23J  Effects  of  Cluster/Carrier-Gas  Collisions 

The  flow  of  gas  that  expands  from  the  supersonic  nozzle  contains  clusters,  carrier 
gas,  and  unclustered  thrust-gas  atoms.  Conceivably,  these  collisicms  could  destroy 
clusters,  alter  their  momentum  distribution,  and  cause  other  effects  that  could  affect  both 
the  quality  and  current  of  the  extracted  clu^  beam.  In  this  section  we  briefly  discuss 
these  processes,  and  we  conclude  that  the  principal  process  of  concern  is  that  of 
momentum-changing  collisions  witMn  the  beam-extraction  system.  We  do  not  at  present 
know  tiie  extent  to  which  these  effects  would  act  in  an  operational  tiuoster,  and  ad^tional 
analytical  study  is  cleaily  needed. 
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In  a  previous  study^  we  found  that  atom-cluster  collisions  are  almost  completely 
elastic  in  character,  and  we  showed  that  the  energy  W(z)  of  clusters  varies  with  axial 
coordinate  z  as 


-{— ) 

W(2)=W,e  ,  (2.3-5) 

where  Wq  is  the  initial  value  of  the  energy.  We  also  showed  that  the  e-folding  distance  zq 
is  given  by 


*•  = 


(2.3-6) 


where  p  is  the  mass  density  of  the  gas  of  mean  atomic  mass  m,  and  Va  is  the  volume  of  an 
atom  (or  mdecule)  of  thrust  gas.  In  evaluating  equation  2.3-6  numerically,  we  find  that  the 
losses  in  energy  can  be  very  substantial:  for  example,  a  cluster  composed  of  300  CO2 
molecules,  passing  through  a  cloud  of  carrier  gas  having  a  (typical)  number  density  of  S  x 
10^  m'3  will  fall  to  1/e  of  its  initial  energy  in  just  40  nun.  llus  energy  is  converted  into 
increased  temperature  of  die  carrier-gas  clmid,  leading  to  a  higher  rate  of  carrier-gas  aimns 
returning  to  collide  with  the  external  side  of  the  decel  electrode  and  the  downstream  side  of 
the  accel  grid  While  these  collisions  will  not  affect  thrust,  they  will  certainly  increase  die 
rate  of  charge-exchange  erosion  of  the  accel  grid  over  the  value  that  would  obtain  without 
the  cluster  collisions.  The  in^xiitance  of  diis  effect  on  accel  erosion  must  be  investigated  in 
future  work. 
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3.0  EXPERIMENTAL  STUDY 


In  this  section  we  describe  the  apparatus  that  we  used  in  our  investigation,  and  we 
presoit  die  expoimental  results. 


3.1  TEST  CONFIGURATION 

Figure  1  was  shown  eariier  to  illustrate  the  basic  arrangement  of  the  experimental 
apparatus;  the  figure  is  repeated  here  for  the  cmivenience  of  the  reader.  The  basic  system 
consists  of  a  plenum,  into  which  thrust  and  carrier  are  premixed;  a  fast  electromagnetic 
valve;  a  supersonic  nozzle;  an  ionizer,  the  beam-extraction  system,  and  a  retarding-potoitial 
analyzer  (^A)  that  is  the  primary  beam  diagnostic.  The  qiparatus  is  basically  the  same  as 
in  Reference  (1),  except  diat  we  have  enlarged  the  ionizer  to  accommodate  the  entire  beam, 
and  we  have  mounted  the  RPA  on  a  movable  arm  to  facilitate  prolong  the  oitire  beam. 


^MANIFOLD 


EXTRACTION 

ELECTRODES 


MOVABLE 

FARADAY 

CUP 


PULSED  VALVE 


OSaLLOSCOPE 


ELECTRON 

EMITTER 


BEAM-EXTRACTION 

ELECTRODES 


Hgure  1,  Duplicated  firom  earlier  in  the  report 

3.1.1  Pulsed  Valve 

In  application  as  an  ionized-cluster  thruster,  continuous-flow  toleration  would  be 
desired  to  provide  maximum  thrust  and  to  minimize  any  losses  that  may  occur  during  the 
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start-up  transients.  However,  continuous  operation  produces  a  gas  flow  rate  of  around 
30,000  liters/s  at  a  pressure  of  100  mPa.  In  space  this  and  larger  flow  rates  would  pose 
no  special  problems,  but  in  a  ground  test  chamber  it  would  require  a  very  large  pumping 
system.  To  avoid  these  difflculdes  we  have  chosen  to  operate  our  system  in  a  pulsed 
r^e,  in  which  the  electrcnnagnetic  valve  is  pulsed  open  for  only  a  few  ms  (typically  0.5  to 
2.0  ms).  Our  pulsed  valve,  which  is  similar  to  one  developed  by  Hagena,^  is  capable  of 
opening  in  as  short  a  time  as  10  |is.  We  determined  in  our  previous  study^  that  the  gas- 
flow  condititms  approach  their  steady-state  values  very  eariy  in  a  gas  pulse. 

3.1.2  Supersonic  Nozzle 

The  supersonic  nozzle  that  we  have  used  has  a  minimum  opening  of  125  pm  and  a 
cone  angle  of  8.5  degrees.  The  cone  is  25-mm  long.  We  found  during  testing  that  the 
skimmer  shown  in  Hg.  1  had  little  or  no  effect,  and  most  of  our  experimental  results  were 
obtained  without  it 

3.1.3  Ionizer 

The  ionizer  consists  of  an  anode  consisting  of  a  short  cylindrical  section  of  wire 
mesh  (16-mm  diameter,  100  mm  downstream  of  the  nozzle — an  angle  of  about  9  degrees, 
to  accept  the  entire  8.5-degree-wide  beam),  outside  of  which  is  the  cathode,  a  hoop-shaped 
filament  which  provides  electrons  thermionically.  No  magnetic  field  is  used. 

In  typical  operation,  we  applied  between  40  and  120  V  between  anode  and 
cathode,  obtaining  discharge  currents  of  between  50  and  130  mA. 

3.1.4  Beam<Extraction  System 

The  beam-extraction  system  that  we  have  used  was  designed  to  be  simple  and  easy 
to  construct  It  consists  of  four  flat  grids  made  from  woven  wire  mesh  having  square 
qtenings  1  mm  on  a  side.  The  grids  each  have  a  transparency  of  about  75%,  and  they  are 
aligned  so  that  the  apertures  are  all  paraxial.  This  construction  differs  substantially  from 
conventional  km-thruster  practice. 

3.1.5  Retarding-Potential  Analyzer  (RPA) 

The  principal  instrumentation  that  we  used  to  analyze  the  ionized  cluster  beam  was  a 
simple  RPA,  consisting  of  a  retarding  grid  and  a  collector,  with  an  aperture  facing  the 
beam.  We  used  this  instrument  to  measure  the  flux,  mass  distribution,  and  profile  of  the 
ionized  cluster  beam. 

3.2  TEST  PROCEDURE 

The  primary  thrust  of  our  experiments  under  this  program  has  been  to  determine  the 
potential  of  ionized  cluster  beams  as  a  mechanism  for  medium-speciflc-impulse  space 
propulsion.  In  this  section  we  identify  these  parameters  and  describe  our  methods  for 
measuring  them. 
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3JL1  Thrurter  Parameters 

The  key  parameters  for  pn^ulsion  are  the  total  power  Pin,  die  thrust  F,  the  specific 
inqnilse  Isp,  dw  prcqiellant-utiliution  efficiency  iim<  and  the  electrical  efficiency  iie>  The 
thrust  is  just 


F  =  4;(mv)  =  ih.v^.  (3.2-1) 

at 

where  m,  is  die  output  mass  flowrate  and  v^|  is  the  cluster  exhaust  velocity.  It  can  also  be 
written  in  terms  of  the  beam  parameters. 


F  =  (3.2-2) 

where  Ib  and  Vb  are  the  beam  current  and  voltage,  respectively.  Met  and  q  are  the  mass  and 
charge  of  the  clusters,  and  y  is  a  correction  factesr  that  accounts  for  cosine-law  thrust  loss  if 
diere  is  any  diveigoice  of  the  beam. 

The  specific  impulse  Isp  is  the  total  thrust-system  impulse  (i.e.,  thrust  x  time) 
divided  the  wei^t  (not  mass)  of  the  [sopellant .  It  can  easily  be  shown  to  equal 

(3.2-3) 

I  1  g 

whoe  g  is  the  acceleration  due  to  gravity  at  the  Earth’s  surface. 

The  propellant  utilization  efficiency  is  just  die  fraction  of  the  input  propellant  mass 
that  ends  iqi  in  die  thrust  beam. 


=  (3.2-4) 

m|  qm, 

where  fr  is  the  fractiem  of  the  total  propellant  mass  that  is  thrust  gas,  and  iid,  ilion,  and 
'Hext  8re  efficieiKues  duit  correspond,  respectively,  to  the  fraction  of  thrust-gas  atoms  diat 
beoxne  clustered,  the  firactiem  of  clusters  that  b«x>me  imiized,  and  the  fraction  of  ionized 
clusters  that  successfully  transit  the  extraction  system  to  provide  thrust  The  remaining 
parameto’,  electrical  efficiency  T|e,  is  ju^  die  beam  power  VbIb  divided  by  die  total  power 
Pin  consumed  by  the  thruster  (including  beam  power,  discharge  power,  nonthrusting 
power  to  the  grids  such  as  accel  current  etc.;  in  a  practical  flight  thrust  system,  the 
efficiency  of  die  power  processor  must  also  be  incorporiued). 

(3-2-5) 

■h 
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The  expressions  given  in  this  section  refer  to  a  beam  that  has  a  single  species  with 
constant  charge-to-mass  ratio.  Appendix  2  shows  that  a  small  but  signiflcant  correction 
reduces  perfonmance  if  certain  types  of  charge-to-mass  distributions  are  present 

3.2.2  Measurement  Methodology 

In  order  to  measure  the  overall  propellant-udlizadon  efficiency  ilin>  we  determined 
the  cluster  current,  the  cluster  mass-to-charge  rado,  and  the  corresponding  flowrate  of  gas 
through  the  supersonic  nozzle. 

In  measuring  the  cluster  current  we  used  our  RPA  as  a  non-energy-discriminating 
Faraday  cup.  We  biased  the  retarding  grid  negadve  (to  prevent  secondary-electron 
emission  fr^  the  collector  plate  which  would  have  appeared  as  false  cluster  current), 
effecdvely  using  the  RPA  as  a  Faraday  cup,  and  measured  the  cluster-current  pulse  over  a 
cross  secdon  of  the  ionized-cluster  beam  emerging  from  the  beam-extracdon  system;  by 
integradng  the  current  density  over  the  beam  cross-secdon,  we  obtained  a  value  for  the  total 
cluster  current 

Another  key  parameter  is  the  cluster  mass-to-charge  rado,  Mci/q.  We  exploited  the 
energy-conservadon  reladcxiship 


(3.2-6) 


to  measure  this  parameter.  Here  Vr  is  the  retarding  potendal  applied  to  the  RPA  grid 
which  just  stops  a  cluster  of  given  mass-to-charge  rado.  The  cluster  velocity  Vd  is 
determined  by  a  dme-of-flight  method,  measuring  the  dme  to  traverse  the  known  distance 
from  the  nozde  to  the  grids  and  dming  the  arrival  delay  from  valve  opening  to  current 
detection  in  the  RPA. 


The  prc^llant-udlizadon  efficiency  tim  is  calculated  frx>m  the  output  beam  current 
and  the  known  total  input  mass  flowrate  ih^of  prc^llant  into  the  nozzle; 


(3.2-7) 


We  can  also  determine  the  product  of  the  ionizadon  and  clusterizadon  efficiencies 
by  measuring  the  total  cluster  current  upstream  of  the  extracdon  system  and  comparing  it  to 
the  input  mass  flow  oi  thrust  gas: 


(3.2-8) 
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3.3  EXPERIMENTAL  RESULTS 


The  technical  issues  outlined  in  preceding  sections  show  that  clusterization, 
ionization,  and  beam  extraction  are  all  key  processes  in  achieving  acceptable  perfoimance  in 
a  thruster.  However,  clusterization  efficiency  and  the  impact  of  carrier  gas  on  propellant- 
utilization  efficiency  are  the  issues  onto  which  we  focused  our  energies,  because  no  amount 
of  improvement  in  ionization  and  beam-extraction  technologies  can  compensate  for  poor 
clusterization  or  a  requirement  for  excessive  carrier  gas. 

We  first  tested  pure  inert  gases  (Ar  and  Xe),  then  inert-gas/carrier-gas  mixtures, 
and  findly  molecular-gas/carrier-gas  mixtures.  We  found  that  pure  inert  gases  cluster  very 
poorly,  inert-gas/carrier-gas  mixtures  are  far  better  but  sdll  not  attractive,  and  molecular- 
gas/carrier-gas  mixtures  show  great  promise.  The  following  sections  describe  these  results 
in  detail 


Hgure  6  illustrates  the  clustering  efficiency  of  the  pure  inert  gases  argon  and  xencm. 
These  data  show  such  low  clusterization  efficiencies  as  to  be  useless  for  practical  thrusting 
systems. 


Figure  6.  Qustering  efficiencies  (ratio  of  clustered  gas  flow  to  total  gas  flow)  measured 
with  single  noble  gases.  The  observed  efficiencies  arc  even  lower  than  the  theoretical 
predictions  (see  Figure  4),  confirming  that  sin^e  gases  are  not  suited  for  cluster  thrusters. 
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Figures  7  and  8  show  the  clusterization  efficiencies  that  we  obtained  with  several 
gas  mixtures,  as  functions  of  the  partial  pressures  of  the  gases.  In  the  data  in  Fig.  8,  we 
have  held  the  stagnation  pressure  Pq  constant  at  825  kPa,  while  varying  the  partial 
pressures  of  the  constituent  gases;  as  thrust-gas  partial  pressure  is  increased,  carrier-gas 
partial  pressure  is  decreased  to  keep  to  total  pressure  the  same  for  each  data  point.  In  the 
upper  half  of  the  figure,  hydrogen  is  used  as  the  carrier  gas,  while  helium  is  used  for  the 
data  shown  in  the  lower  half.  Remembering  the  arguments  associated  with  Fig.  5,  we  see 
that  the  lower  clusterization  fractions  and  less  desirable  masses  of  inert-gas  clusters  derived 
either  from  pure  gases  or  from  thrust-gas/carrier-gas  mixtures  drive  us  in  the  direction  of 
the  nx>lecular-gas/carrier-gas  mixtures.  Accordingly,  we  devoted  the  remainder  of  our 
experimental  effort  to  these  promising  mixtures. 
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Figure  7.  Qustering  efficiencies  measured  with  gas  mixtures.  At  pressures  near  27  kPa  some 
cluster  fonning  gases  are  seen  to  be  clustered  nearly  completely.  Accmdingly,  these  gases  are  well 
suited  for  use  in  dmistos.  However,  diey  remain  ^ficient  in  cme  respect:  They  need  a  larger  than 
»pectBd  quantity  of  carrier  gas  (close  to  800  kPa  with  27  kPa  thrust  gas).  This  leads  to  a  lower 
duin  predicted  propellant  utilization.  The  highest  measured  utilizsaion  value  readied  42% . 

While  the  data  in  Fig.  7  show  die  effoct  of  thrust-gas  partial-pressure  variation  of  fixed  tool 
pressure.  Fig.  8  illustrates  d^  effects  of  variations  in  total  pressure,  with  thrust-gas  fraction  as  a 
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parameter.  It  is  clear  that  higher  total  (=stagnation)  pressures  produce  superior  clusterization, 
except  for  the  cases  in  which  higher  thrust-gas  partial  pressures  are  used.  The  best  clusterization 
show  in  the  figure  is  about  95%. 

Two  imptmant  results  are  evident  in  Rgure  7;  first,  the  carbon-dioxide  and  ethylene 
(C2H4)  data  show  clusterization  efficiencies  approaching  100%,  and  second,  there  appears 
to  be  an  optimum  ratio  of  thrust  gas  to  carrier  gas  at  this  value  of  the  stagnation  pressure. 

An  obvious  speculation  that  could  explain  this  effect  is  that  thrust-gas-atom  collisions  with 
clusters  (the  collisions  that  cause  clusters  to  grow)  are  too  infrequent  under  conditions  of 
low  thrust-gas  partial  pressure,  and  the  carrier-gas  fraction  is  insufficient  to  cool  the 
growing  clusters  under  conditions  of  high  thrust-gas  partial  pressure. 
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Figure  8.  Dependence  of  clustering  efficiencies  on  the  total  (stagnation)  pressure,  widi  the 
ratio  of  clustering-to-canier  gas  pressures  as  a  parameter. 

From  the  standpoint  of  thruster  performance,  our  best  data  were  obtained  with  the 
mixture  C02:H2  =  27 1^8:800  kPa.  Table  1  illustrates  the  performance  calculated  from  the 
measurables  identified  earlier  in  this  section.  There  are  several  important  conclusions  that 
can  be  reached  studying  this  table.  The  fraction  of  the  thrust  gas  that  actually  clustered 
is  quite  high,  94%,  but  tluit  fraction  represents  only  42%  of  the  total  mass  input  (the 
reniainder  being  carrio’  gas).  Also,  the  prc^Uant  utilization  efficiency  is  only  15.0%,  the 
reduction  from  42%  occurring  in  the  woven-wire-mesh  beam-extraction  system  (deduced 
fnom  measurements  of  the  ionized-cluster  current  on  both  sides  of  the  extraction  system). 
This  low  value  of  propellant  udlizadon  efficiency  brings  the  specific  impulse  to  203  s 
(uncorrected  for  beam  divergence  and  the  fxesence  of  a  disoibudon  in  q/m;  if  qm  had  been 
100%  and  there  were  negligible  beam  divo-gence  and  only  a  single  qAn  q)ecies,  Isp  would 
have  been  1360  s).  These  observadons  point  out  that  getting  the  carrier-gas  flow  down 
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substantially  will  be  necessary  before  ionized-cluster  beams  have  a  chance  in  the  space- 
pix^ulsicm  arena.  The  origin  of  the  severe  current  losses  incuned  in  passing  through  the 
extraction  system  are  unknown,  but  it  could  reside  in  the  primitive  nature  of  our  laboratory 
extraction  system  (or  possibly  the  presence  of  the  retarding  grid),  or  it  may  be  only 
i^jparent,  due  to  charge  exchange  between  the  cluster  ions  and  the  large  amount  of  neutr^ 
carrier  gas.  This  questicm  requires  further  investigation. 


Table  1.  lonized-Cluster  Thruster  Performance 

C02:  27  kPt,  H2:  800  kPt 
7.4x1 0-«  kg/t 


Gat  mixture 
Total  matt  flow 
Thruat  maai  flow 
Cluiter  Current 
Chiater  formation  and  ionization  efficiency 
Propellant  utilization 
Average  cluaier  tize 
Extracted  chiater  current 
Beam  voltage 
Ionization  voltage 
Ionization  cntrettt 
Berm  power 
Ionizer  power 
Electrical  efficietKy* 

Thrust 

ThruaUPower  ratio 
Speciflc  impulae/y 


3.3x10-*  kg/s 

23.2  mA 

94  %  (42%  of  total  gas  mass) 

IS.0%  (after  extraction)* 

296  molecules/e 

8.2  mA 
12  kV 
100  V 
100  mA 
98.4  W 
SOW 

663%  (would  be  90.8  %  without  filament) 

14.8  mN 

99.8  mN/kW 
203  s* 


*  illustrates  severe  extraction-system  losses 


3.4  SUPPORTING  MEASUREMENTS 


In  this  section  we  describe  a  number  of  investigations  that  we  carried  out  to  better 
define  our  measurements  of  nozzle-exhaust  characteristics,  the  cluster-ionization  process, 
and  the  cluster  beam. 


3.2.1  Nozzle  Characteristics 

In  Fig.  9(a),  we  show  the  velocity  of  pure  gases  expanding  through  the  supersonic 
nozzle,  as  detennined  by  a  time-of-flight  method.  We  believe  that  the  decrease  in  velocity 
at  low  stagnation  pressures  may  be  due  to  boundary-layo-  effects.  At  higher  stagnation 
pressures,  the  flow  velocities  approach  their  thetmtical  values,  based  tm  the  theory  of 
isentrqiic  expansion  into  vacuum.  In  Fig.  9(b),  however,  we  see  a  more  complex  shqie 
(inset)  of  gas  pressure  arriving  at  the  detector.  We  believe  that  the  three  portions  of  the 
arriving  pulse  cmrespond  to  carrier  gas,  uncondensed  thrust  gas,  and  finally,  clusters.  We 
cannot  tell  from  these  data  whether  the  three-velocity  characteristic  ctmtinues,  or  is 
associated  only  with  the  leading  edge  of  the  pulse. 


22 


3 


PURE  GASES 


Figure  9.  Exhaust  velocities  of  single  gases  (top)  and  of  a  gas  mixture  (bottom),  derived 
from  leading-pulse-edge  time-of-flight  measurements.  In  the  gas  mixture  with  its  high 
cluster  content  three  different  velocity  groups  can  be  distinguished:  fast  H2,  a  slower 
mixture  of  Ar  and  H2,  and  still  slower  Ar  clusters.  It  is  uncertain  whether  these  velocities 
are  just  leading-edge  effects  or  are  steady-state  phenomena. 

In  Fig.  10,  we  contrast  the  predictions  of  isentropic-expansion  theory  with 
measured  data  of  the  velocity  of  a  pure  gas  (helium)  and  an  argon-helium  mixture. 
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STAGNATION  PRESSURE,  kPa 

Figure  10.  Gdculated  and  measured  nozzle  flow  rates.  The  calculadcms  were  based  mi 
isentropic  flow  conditions  in  die  nozzle  throat 

3.2.2  Ionizer  Characteristics 

Figure  1 1  illustrates  the  characteristics  of  our  simple  ionizer.  As  ionization  current 
is  increased,  die  cluster  size  is  unaffected,  suggesting  that  the  ionizer  does  not  cause  cluster 
fragmentation;  and  the  imiized-cluster  current  rises  and  saturates,  suggesting  that  we  are 
reaching  eidier  complete  ionization  of  all  clusters,  or  that  some  odier  saturation  mechanism 
is  coming  to  play. 

In  Fig.  12,  we  show  the  effects  of  varying  the  bias  mi  the  grid  nearest  the  plasma 
(which  is  intended  to  prevent  single  ions  from  entering  the  beam).  The  sharp  comer  in  the 
characteristic  corresponds  to  the  ionizer  anode  potential  of  100  V,  illustrating  that  the 
potential  cf  the  plasma  near  this  grid  resides  very  close  to  the  ionizer  anode  potential. 

As  discussed  in  earlier  sections,  our  simple  ionizer  requires  excessive  power  in  the 
form  of  filament  heat,  and  has  other  characteristics  that  produce  non-<^timal  perfmmance. 
We  feel  that  a  hollow-cathode-based  Penning-type  discharge  offers  Ae  best  promise  fm* 
averting  these  difficulties,  but  program  resources  did  not  permit  us  to  investigate  diis 
possibility.  We  did,  however,  test  a  cold-cathode-based  ionizer,  which  we  found  to  be 
even  less  effidoit  than  our  filament  ionizer. 

The  cold-cathode  ionizer  consisted  of  a  cylindrical  cone  surrounding  the  expanding 
cliater  beam,  inliich  was  biased  to  between  -500  V  and  -3000  V  relative  to  the  first  (single¬ 
ion-retarding)  grid  of  the  beam-extraction  system.  Figure  13  compares  the  cold-cathode 
kmizer  perfmmance  with  that  of  the  filammtt-based  ionizer. 
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3.2.3  Beam  Extraction  System  Characteristics 

As  discussed  earlier,  we  incurred  a  substantial  perfannance  penalty  because  of  the 
characteristics  of  the  beam-extraction  system.  Late  in  the  program,  we  investigated  the 
possibili^  that  the  problem  may  have  been  that  die  beam  was  reaching  space-charge  limits 
within  the  extraction  system.  This  may  have  been  part,  but  is  certainly  not  all,  of  the 
problem. 

For  our  eariy  data,  we  used  grid-to-grid  gaps  of  6  to  8  mm  in  the  extraction  system. 
To  determine  whether  these  large  gaps  were  limiting  the  extracted  cluster  currents,  we 
reduced  the  spacings  to  2.5  mm,  and  found  the  effects  shown  in  Fig.  14.  Since  the  space- 
charge-limited  current  between  plane-parallel  electrodes  goes  as  the  inverse  of  the  square  of 
the  interelectrode  gap,  theory  would  predict  an  improvement  of  a  factor  of  6  to  8,  rather 
than  the  modest  factor  achieved. 

Figure  15  is  a  radial  proHle  of  the  current  density  emanating  from  the  beam- 
extraction  system.  The  domed  profile  is  commonly  seen  in  ion  thrusters,  where  the  highest 
plasma-generation  rate  is  near  the  center,  but  no  comparable  mechanism  is  likely  to  be 
present  in  our  ionizer.  This  effect  and  the  anomalous  cluster  loss  through  the  extraction 
system  may  be  due  to  aperture-to-aperture  alignments  that  are  inappropriate  for  heavy 
clusters. 


IONIZER  CURRENT,  mA 


Figure  1 1.  Quster  size  and  current  as  a  function  of  ionizer  current  The  cluster  current 
iqrproaches  a  limit  value  which  is  taken  to  be  full  ionization.  The  independence  of  cluster 
size  on  ionizo’  current  is  taken  as  an  indication  that  multiple  ionizations  do  not  increase  witii 
current 
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Figure  12.  Retarding  potential  characteristic  of  the  cluster  current,  measured  at  the 
retarding  grid  next  to  the  ionizer.  The  sharp  break  of  the  curve  at  anode  potential  indicates 
that  the  potential  within  the  ionizer  is  uniform,  and  that  the  plasma  potential  is  very  close  to 
aiKxle  potential 
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Figure  14.  Observation  of  space-charge-limited  flow  in  the  beam-extraction  system. 
Initially,  the  intergrid  gap  was  8  mm  wide  and  could  not,  even  at  the  highest  voltages, 
transmit  the  available  cluster  currents.  The  situation  improved  significantly  after  the  gap 
width  was  decreased  to  2.5  mm.  The  dashed  curve  shows  the  substantial  magnitude  of 
secondary  electron  currents.  These  were  eliminated  at  the  Faraday  cup  with  the  help  of  a 
slightly  negatively  biased  grid. 
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Figure  15.  Transverse  cluster  beam  distribution.  Beam  profiles  were  integrated  to 
determine  die  total  beam  current 
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4.0  SUMMARY  AND  CONCLUSIONS 


Under  this  program,  we  have  taken  the  first  steps  in  demonstrating  that  space 
propulsion  using  ionized  clusters  is  feasible,  although  a  great  deal  of  scientiHc  and 
engineering  work  stands  between  our  present  status  and  having  available  flightready 
cluster-ion  thrusters.  We  have  demonstrated  that  we  can  form  beams  of  clustered  atoms  or 
molecules,  ionize  dtem  and  extract  them  into  a  beam,  achieving 

•  Thrust  levels  >10  mN  at  practical  beam  voltages 

•  Qusteiizadon  efficiencies  >90%  with  desirable  thrust  gases 

•  Ionization  efficiencies  >90% 

Despite  these  accomplishments,  we  are  presently  hampered  by  several  performance- 
robbing  processes: 

•  Excessive  amounts  of  earner  gas  are  needed  for  high  clusterization,  degrading 
mass-utilization  efficiency  to  40%. 

•  There  is  a  substantial  drop  in  cluster-ion  current  after  passing  through  the 
beam-extraction  system;  this  further  degrades  mass-utilization  efficiency  to 
(Hily  15%. 

•  The  ionizer  requires  excessive  power  (50  W),  reducing  electrical  ^ciency 

We  feel  that  these  negatives  are  challenging  but  pose  no  clear  baniCTS  to  the  future 
development  of  cluster-ion  propulsion  to  a  practical  r^ty.  Without  further  study,  we 
cannot  with  certainty  fmesee  the  path  that  will  lead  to  a  resolutiem  of  these  issues,  but  we 
outline  some  ideas  below  that  we  believe  wmthy  of  investigatitm.  Some  of  these  ideas 
require  scientific  study,  and  odiers  are  simply  matters  of  engineering  that  we  include  for 
completraess. 

•  Carrier-Gas  Injection.  The  carrier  gas  is  presently  injected  through  the 
supersemk  nozzle,  as  a  mixture  with  the  thrust  gas,  altiiough  the  cooling 
benefits  of  the  carrier  gas  are  not  needed  until  after  the  saturation  point  has 
been  reached.  After  reviewing  our  present  data,  we  believe  that  it  may  be 
boieficial  (i.e.,  require  less  carrier  gas)  to  introduce  the  carrier  gas  stmiewhat 
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downstream  of  the  nozzle  orifice  so  that  it  does  not  undergo  the  supersonic- 
flow  transformation. 

•  Carrier-Gas  Flow  Control.  Presently  the  drift  region  between  the  nozzle 
and  the  ionizer  is  completely  open  to  the  vacuum  chamber,  a  design  that  is 
intended  to  avoid  Paschen  brealcdown  in  the  beam-extraction  system.  This 
design  achieves  its  objective,  but  it  may  be  possible  to  substantially  enclose  the 
drift  region  without  incurring  breakdown.  This  measure  will  inhibit  the  early 
loss  of  carrier  gas  and  permit  some  of  it  to  be  recycled  from  the  enclosure 
walls.  This  may  be  more  important  in  a  continuously  operating  cluster-ion 
thruster  than  it  is  in  our  pulsed  device. 

•  Magnetic  Filter  to  Replace  the  Retarding  Grid.  The  retarding  grid  is 
an  inefficient  design  for  several  reasons:  it  collects  electrons  from  the 
discharge,  causing  an  unnecessary  power  drain,  it  probably  distorts 
trajectories  of  cluster  ions  entering  the  beam-extraction  system,  and  it  almost 
certainly  interferes  with  the  formation  of  the  necessary  meniscus-shaped 
plasma  sheath  that  forms  over  each  screen-electrode  extraction  aperture  in 
beam-extraction  systems  that  lack  retarding  grids.  Magnetic  filters,  composed 
of  strong  magnetic  fields  transverse  to  the  flow  direction,  have  been 
demonstrated  to  effectively  separate  species  with  different  charge-to-mass 
ratios,  a  capacity  that  electrostatic  filters  normally  lack  (only  the  velocity 
difference  Iwtween  the  cluster  ions  and  light  ions  allows  them  to  be  separated 
by  electrostatic  means).  Such  a  filter  could  replace  the  retarding  grid, 
improving  the  performance  of  the  extraction  system  and  avoiding  the  power 
drain. 

•  Magnetic  Confinement  of  the  Discharge  Plasma.  The  discharge 
plasma  is  presently  not  confined  at  all,  causing  what  are  probably  undesirably 
low  plasma  densities  at  the  extraction  system,  and  wasting  most  of  the  energy 
of  the  100-V  ionizing  electrons.  By  using  boundary  magnetic  fields  to  confine 
the  plasma,  these  problems  could  directly  avoid^,  and  the  added  discharge 
efficiency  could  conceivably  allow  operation  at  voltages  low  enough  that 
multiply-charged  clusters  (with  their  associated  performance  penalty)  are  not 
formed. 

•  Hollow  Cathode  to  Replace  the  Filament.  This  obvious  measure  will 
save  a  large  amount  of  the  power  that  is  presently  used  to  heat  the  filament 
emitter  (around  3  W  vs.  40  W).  Hollow  cathodes  also  have  the  lifetime 
required  of  flight  thrusters. 
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APPENDIX  1.  SYMBOLS  USED  IN  THIS  REPORT 


lip  Specific  Impulse 

Mp  Total  Pnqiellam  Mats 

Total  Efficiency 

y  Beam-Divergence  Comcdon  Factor 

g  Acceleration  due  to  Earth’s  graviq^ 

F  Thrust 

P  Pdwer 

Ir  Total  Mission  Impulse  Requirement 

V  Ratio  of  SpeciEc  Heals 


'P 

Cy 

T 

V 

Vd 

L 

n 

j* 

6o 

ZQ 

Vd 

m, 

e 

Md 

q 

Vg 

w. 

m 

Vb 

Wd 

Wo 

p 

V. 

*# 


Specific  Heat  at  Constant  Pressure 
Specific  Heat  at  Constant  Volume 
Absolute  Temperature 
Vohnne 

Cluster  Velocity 
Lengdi  of  ionizer  cell 
Ionization  Cross-Section 

Oas  Number  Densi^ 

Electron  Current  Density  In  Ionizer 

Permitliviqr  Of  Free  Space 

Orid-To-Grid  Gqp 

Disdiarge  Voltage 

Maas  Of  The  Electron 

Charge  Of  The  Electron 

Total  Cluster  Mass 

Charge 

Voltage  On  Retarding  Grid 

Energy  Of  Atom 

Mass  Of  Atom  Or  Ion 

Maximum  Adiievdtle  Specific  Impulse 

Beam  Voltage 

Cluster  Energy 

Initial  Chisier  Energy 

Gas  Mass  Density 

Atomk  Volume 

Rate  Of  OMput  or  Mass  fa  Thrust  Beam 


h 

n, 

r\m 

h 

lUd 

lliaa 

R 

N 

PW 


Rate  Of  bpiu  Of  Mass  To  Thfuster 

BeamCuimt 
Electrical  EfBcieiicy 
Propellant  Mau-Utilizatkm  Efficiency 
Man  Fraction  Of  Thrutt  Gu 
Efficiaicy  Of  Cluttering  Tfanut  Gas 
Efficiency  Of  Ionizing  Clusten 
Beam  Voltage/Total-Extraction  Voltage 
Number  Of  Mdecukt  In  Chisier 
Total  bqnt  Pbwer  To  Tfanitter 
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APPENDIX  2.  EFFECT  OF  MASS-TO>CHARGE  DISTRIBUTION 
ON  THRUSTER  PERFORMANCE 


As  mentioned  in  Section  2.3,  the  presence  of  a  distribution  of  charge-to-mass  ratios 
q/m  in  the  beam  can  cause  a  thruster  performance  loss.  This  loss  occurs  because  specific 
impulse  and  the  thrust-to-power  ratio  are  related  in  a  way  that  depends  upon  the  distribution 
of  q/m  values.  In  this  appendix  we  derive  expressions  for  specific  impulse  and  thrust-to- 
power  ratio  that  are  applicable  to  situations  in  which  a  distribution  of  q/m  values  is  present, 
and  we  show  the  impact  of  the  presence  of  the  distribution. 

To  summarize  our  results,  we  find  that  the  specific  impulse  of  a  beam  with 
distributed  qAn  is  the  same  that  would  obtain  if  the  beam  consisted  of  uniform  q/m  at  the 
mean  value  of  the  distribution.  For  a  given  specific  impulse,  the  thrust  to  power  ratio  is 
reduced  by  a  fraction  that  is  of  order  o?/3ti^,  where  o  is  the  standard  deviation  of  the 
distribution,  and  q  is  its  mean. 

A2,l.  SPECIFIC  IMPULSE 
Specific  impulse  is  defined  as 


I,  =  FAt/Mg,  (A2-1) 

where  F  is  tiie  thrust.  At  is  the  duration  of  thrusting,  g  is  the  acceleration  due  to  the  Earth’s 
gravity,  and  M  is  the  total  propellant  mass.  Recognizing  M/At  =  ih  /  q.,  where  ih  is  the 
mass  flowrate  in  the  thrust  beam,  and  qm  is  tiie  propellant-utilization  efficiency,  we  have 


The  thrust  F  is  given  as 


F  =  ihv  =  niY 


(A2-2) 


(A2-3) 


and  we  define 


S  .S 


(A2-4) 

Vm 

giving 

F  =  (A2-5) 

If,  however,  the  mass  flowrate  is  distributed  over  many  ^  with  a  distribution 
function  f(^,  then 


dF  =  iiif(5)^T5«l6, 
and,  assuming  y  to  be  independent  of  4, 

F  =  mV2STp(4«l5- 

Finally,  using  A2-2, 

9  9 

In  equation  A2-8,  note  that  the  integral  is  sinq>ly  the  average  value  (^)  =  ^ 

that  the  usual  single-q/m  equation  holds  if  one  substitutes  an  appropriately  weighted 
average  value  of  q/m. 

A2J  THRUST-TO-POWER  RATIO 

We  have  already  obtained  an  expression  for  the  thrust,  eqn.  A2-7,  so  now  we  need 
the  power.  In  the  single-qAn  case,  the  beam  power  is  simply  VbIb>  and  Is  is  given  by 

I  =  (A2-9) 

m 

die  case  where  ^  is  distributed,  the  beam  power  is  still  VbIb;  however,  the  beam  current 
now given 


(A2-6) 


(A2-7) 


(A2-8) 


I,  =  *J^5’fft)d5.  (A2-10) 

Coadnning  A2-10  and  A2-7,  we  have 
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F/p=ri. 


\v. 


(A2-11) 


Again  notice  that  the  numerator  is  the  average  (first  moment)  of  the  distribution 
f(4),  and  denominator  is  its  second  moment  Equation  A2-1 1  can  therefore  be  rewritten 
more  conq)actly. 


F/P  = 


(A2-12) 


The  first  moment  is  the  mean  s  but  the  second  moment  differs  from  the 

second  central  moment  or  variance,  o^: 


(A2-13) 

By  expanding  the  right-hand  side  of  Equation  A2-13  and  rearranging  the  terms,  we 

obtain 


{5’)  =  o'+V.  (A2-I4) 

from  which  can  expand  the  fraction  -^r  in  equation  A2-12  (assuming  0  «  q)  to 

obtain 

The  first  Gowest-order)  term  in  equation  A2-15  resembles  the  thrust-to-power 
expression  fm*  a  thruster  involving  a  single  charge-to-mass-ratio  species;  the  subsequent 
terms  give  ctnrections  when  there  is  a  spread  in  the  charge-to  mass  ratios. 

In  this  program  we  were  unable  to  accurately  measure  the  distribution  of  charge-to- 
mass  ratios,  so  we  presently  have  no  numbers  to  put  into  equation  A2-1S  to  estimate  our 
performance  loss.  However,  one  thing  is  clear  from  Eqn.  A2-1S:  if  the  standard  deviation 
o  of  the  charge-to-mass  distribution  were,  say,  10%  of  the  average,  ii,  the  thrust  loss  due 
to  q^qnead  will  be  less  than  1%.  This  situation  might  arise  if,  for  example,  all  clusters 
had  the  same  charge  and  a  10%  variation  in  mass.  A  potential  problem  arises  when  there 
are  multiply-charged  clusters.  Suppose  that  the  cluster  beam  consisted  of  equal  parts  of 
singly-ionized  and  doubly-ionized  clusters  having  mass  250.  In  this  hypothetic  exan:q)le, 
we  find,  t]  =  6x10*3  and  o  =  1.41x10*3,  in  units  of  e/Mii,  where  Mh  is  the  mass  of  the 


36 


hydrogen  atom.  The  loss  in  thrust-to-power  ratio,  to  third  order,  is  then 
o^/n^  -  5.3  X  10-2,  or  5.3%. 

For  readers  with  ion-propulsion  backgrounds,  the  foregoing  results  may  seem 
strange,  since,  for  example,  the  presence  of  doubly  charged  ions  in  an  ion-thruster  beam  is 
normally  thought  to  cause  a  substantial  performance  loss.  The  resolution  of  this  apparent 
contradiction  is  that  the  analysis  given  in  this  appendix  has  determined  the  performance  loss 
due  to  the  fact  that  the  variance  of  the  distributicm  of  charge-to-mass  ratios  in  the  propellant 
is  non-zero,  and  effects  of  variation  of  the  mean  of  the  propellant  mass  are  not  discussed. 
In  conventional  ion-propulsion  practice,  thruster  parameters  are  calculated  with  the 
assumption  that  ail  of  the  propellant  is  of  a  single  charge-to  mass  ratio  (e.g.,  Xe'*'),  even 
with  the  foreknowledge  that  other  species  will  be  present  (e.g.,  Xe'*'*').  After  the  fraction 
of  Xe'*"*’  is  determined,  the  performance  that  was  calculated  assuming  that  the  beam  was 
(xxnposed  entirely  of  Xe***  is  conected  to  reflect  the  tme  beam  composition.  This  correction 
factor  appears  to  be  a  performance  loss  (decrease  in  thrust-to-power,  increase  in  specific 
impulse),  but  in  fact  it  is  simply  the  effect  of  using  a  propellant  of  higher  average  charge-to- 
mass  ratio  than  pure  Xe'*';  it  is  similar  to  the  effect  of  substituting  Kr  (a  lighter  gas)  for  Xe 
as  a  propellant,  and  has  nothing  to  do  with  the  fact  that  multiple  species  per  se  are  present.. 
(If  the  tlmister  had  been  assumed  to  operating  on  pure  Xe**"**  instead  of  pure  Xe***,  then  the 
correction  factor  accounting  for  the  presence  of  Xe'*'  would  appear  to  be  an  increase  in 
thrust-to-power  ratio.) 

Since  the  equations  governing  specific  impulse  and  thrust-to-power  ratio  contain  the 
charge-to-mass  ratio  only  as  the  product  qVe/m,  any  change  in  q/m  can  be  compensated  by 
altering  Vb  to  keep  the  product  constant.  A  thruster  with  Tim  =  T)e  =  1  will  then  exhibit 
identical  performance  with  higher  q/m  and  lower  Vb-  A  real  ion  thruster  may,  of  course, 
incur  a  change  in  electrical  or  propellant-utilization  efficiency  when  the  beam  voltage  is 
changed,  so  that  the  overall  performace  may  indeed  change  as  a  result  of  the  correction. 
This  change  is  a  different  topic  entirely,  having  to  do  with  thruster  performance  curves 
rather  than  the  effect  of  distributed  charge-to-mass  ratios. 
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